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Abstract 

The authors of [RB| have observed the following remarkable phe- 
nomenon during their experiments. If two oppositely charged droplets 
of fluid are close enough, at first they attract each other and touch even- 
tually. Surprisingly after that the droplets are repelled from each other, if 
the initial strength of the charges is high enough. Otherwise they coalesce 
and form a big drop, as one might expect. 

We present a theoretical model for these observations, using mean 
curvature flow. The local asymptotic shape of the touching fluid droplets 
is that of a double cone, where the angle corresponds to the strength of 
the initial charges. Our model yields a critical angle for the behaviour of 
the touching droplets and numerical estimates of this angle agree with the 
experiments. This shows, contrary to general belief (see |RB] and |RB1) 1. 
that decreasing surface energy can explain the phenomenon. 

To determine the critical angle within our model we construct appro- 
priate barriers for the mean curvature flow. In |AI) Angenent, Chopp and 
Ilmanen manage to show the existence of one-sheeted and two-sheeted, 
self-expanding solutions with a sufficiently steep double cone as an initial 
condition. Furthermore they provide arguments for nonuniqueness even 
among the one-sheeted solutions. We present a proof for this, yielding a 
slightly stronger result. Using the one-sheeted self-expanders as barriers 
we can determine the critical angle for our model. 

1 Introduction 



In |RB| the behaviour of oppositely charged droplets of fluid is investigated. 
Droplet motion induced by electrical charges occurs in a vast number of applica- 
tions, including storm cloud formation, commercial ink-jet printing, petroleum 
and vegetable oil dehydration, electrospray ionization for use in mass spectrom- 
etry, electrowetting and lab-on-a-chip manipulations (see also !RB ). 

The phenomenon can be described as follows. Two close enough oppositely 
charged droplets of fluid attract each other and converge. Both, experiments and 
numerical simulations (see [RBI] ), provide evidence that a short-lived bridge 
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is formed between the droplets, which instantly causes the charges to be ex- 
changed. The bridge between the touching droplets has the local asymptotic 
shape of a double cone. Furthermore the charges determine the angle of the 
double cone, where a lower charge corresponds to a steeper cone (i.e. having a 
larger acute angle with the rotation axis). 

We want to study the behaviour of the system after the droplets have 
touched. One might think that coalescence occurs and that one big drop is 
formed. However experiments in |RB| have shown that above a critical field 
strength the droplets do not coalesce after touching but are repelled from each 
other. 

Following an idea by P. Topping, we present a theoretical model for this phe- 
nomenon. It is assumed that after the two droplets have touched and exchanged 
their charges, the motion of the system is driven by minimization of energy. To 
model this mathematically we use the mean curvature flow which is the gradient 
flow of the surface area (see e.g. |EC] V Our model has two main advantages 
over the theoretical approaches presented in RB j and [RBI] . Firstly we are 
able to show that minimization of surface energy can explain the observations 
from the experiments. Secondly our results are mostly independent from any 
assumption on the precise local conical shape formed by the touching droplets. 

We define the double cone in R™ for < a < ^ as 

D a = jx £ M" : |(x 2 , . . . , x n )\ 2 — tan 2 a ■ x\~^ . 

The results we obtain from our model can be summarized as follows. 

Assume that two initially oppositely charged droplets of fluid after touching 
have the local shape of a smoothing of a double cone D a in K 3 . Also assume 
that their motion is governed by minimization of area which we model using the 
mean curvature flow. 

Then there is a critical angle a CT n with the following properties. If a < a cr u 
the droplets are repelled from each other. If the associated smoothing lies outside 
of D a and a > a cr n the droplets coalesce and form one big drop. 

Using appropriate barriers and a level-set flow argument we can conclude 
that a cr u is precisely the critical angle for the existence of one-sheeted, self- 
expanding evolutions of D a C K 3 , which means a cr u ~ 66°. 

These formulations are made precise in the sequel. The critical angle of 
60° — 70°, observed during experiments (see |RBj ). agrees with our prediction. 

A family of smooth, immersed hypersurfaces (M t ) teI (I a real interval) in 
M. n is called a solution of the mean curvature flow if 

3x 

— =H, xeM t ,tel. (l) 

H = —Hv is the mean curvature vector and v a choice of unit normal. ([T]) is 
equivalent to ^ = Am,i (Ajvf, is the Laplacian of M t ). 

Based on the ideas in jAIj we present a new proof for the existence of one- 
sheeted self-expanders with the double cone D a for a large enough as an initial 
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condition. We call a solution of ([T]) self-expanding if 

M t = Vt-Mx, te(0,oo). 



(2) 



The singular initial condition D a is here understood to be attained locally in 
the sense of Hausdorff distance. 

Theorem 1. For n > 3 there exists a critical angle a* rit (n) £ (Oj §) with the 
following properties. 

For any angle a > a* rit there exist at least three distinct, smooth, rotationally 
symmetric evolutions of the double cone D a which are self- expanding. Two of 
these evolutions are one-sheeted and one is two-sheeted. 

For a = a* rit at least one one-sheeted and one two-sheeted self- expanding, 
smooth, rotationally symmetric evolution of D a exist. 

Here "one-sheeted" and "two-sheeted" refer to the number of connected 
components of the solutions. The existence of one-sheeted self-expanders was 
first proved in |AI| . Additionally our proof provides nonuniqueness among one- 
sheeted solutions, which was also first stated in |AIj . 

Nonuniqueness here corresponds to fattening of the level-set flow. For a < 
a* rit the two-sheeted evolution of D a is unique and therefore we have non- 
fattening of the level-set flow (see |Alj). 

Using the one-sheeted self-expanders we can show that a cr a = a* rit (3). 

2 Existence of self-expanders 

In this section we want to study the evolution by mean curvature of the double 
cone D a , following [Alj . We are interested in solutions (M t ) t£ ^ , that satisfy 

M t is rotationally symmetric and self-expanding. (3) 

From ([1} and @ one can see that for a solution of the mean curvature flow to 
move self-expanding is equivalent to the self-expanding equation 

H = x £ Ml (4) 

Under condition ([3]) this equation becomes an ODE. 

We write x = (xi, X2, ■ ■ ■ , x n ) = (xi,x) for x € M™ = K x K"" 1 . For a curve 
7 in R 2 , that is symmetric with respect to u i— > — u we define the corresponding 
surface of rotation 

if(7) = {(ii,i)€R" : {x 1} \x\) = (y,u) G 7 }. 

For — 7r < a < tt let a a be the closed ray {i(cos a, sin a) : t > 0} in M 2 . The cone 
in E n with angle a is defined as C a — M (a a U cr_ Q ). We get D a = C a U C T _ Q . 

Following |AIj we get an equation for 7 from ([3]). Parametrizing 7 by ar- 
clength s we define 8 £ [0, 2tt) along 7 by setting 7s = (cos 9, sin 9) for the 
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tangent vector j a . The left-handed unit normal is given by v = (— sin 8, cos 6). 
Let k be the curvature vector of 7 and k — k-T7. Equation <j3j) for M (7) becomes 

/T- ~~ 2 y ii 

fc cos6> + — sin 6* — — cos 8 = 0, (y,u) e 7. (5) 

A solution of this equation creates a smooth surface M (7). 

The following lemma is due to Angenent, Chopp and Ilmanen (see | Alj ) . 
The first part is a consequence of results from Ecker and Huisken on graphical 
mean curvature flow (see |EH] ) . 

Lemma 2.1. (i) (Two-sheeted case) 

For a £ (0,7r) there exists a unique, smooth, connected curve 7(a), solving 
0) and asymptotic to C a . Furthermore, unless 7 is the u-axis, 7 is the graph 
of a positive, convex (or negative, concave) even function y = y(u). 

(ii) (One-sheeted case) 

Let 7 be another smooth, connected curve solving |5[) which meets {u > 0}. 
Then 7 lies in {u > 0} and is asymptotic to a a U a/3, where < a < j3 < it. 
If 7 meets the u-axis at a right angle, then 7 is the graph of a positive, even 
function u — u (y) which is monotone for y and (3 = ir — a. 

We define 

A = |a e ^0, — ^ : 3 connected 7 solving ([3]), asymptotic to a a U ov- Q j 

and a* rit (n) = inf A. 

Now we want to focus on smooth, connected curves 7 solving ([5]) which 
are graphs of even functions (part (ii) of lemma 12. ip . After imposing initial 
conditions we get the following initial value problem for u 

Uyy = (l + (u y ) 2 )(^(u-u y y)+ 1 ^-^y u(Q) = C>0,u y (0) = 0. (6) 

For symmetry reasons it is enough to consider it | [0,00) which we denote again 
by u. For y > we denote by a(y) = arctan^^j the signed angle, that u 

makes with the positive y-axis. First we note the following basic observations 
regarding the solution u. 

i. Every critical point of it is a strict local minimum. 

To see this compute u yy = f + > from ([6]) whenever u y = 0. 



ii. u y > for all y > and every critical point of u y is a strict local maximum 

c I n- 

2 c 

,0) („, \2 



The first part is a direct consequence of [1] and the fact that u yy (0) = § + > 



0. For the second part compute u y ' = — (l + (it y ) 2 ) (n — 2)^f < from @ 
using the first part whenever u yy — 0. 

iii. Every critical point of a (y) is a strict local minimum. Therefore u is asymptotic 
to a ray a a , a = lmij / _ 1 . 00 a (y). 

Let = a v = ^ip- (yu v - u) , y > 0. Using © we get = yu y - u = 
2£n-22 and therefore > But then _ JI^mv > q 
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iv. < a = limy^oo a (y) < ■§. 

This is poved in [AI] using the clearing out lemma (see appendix |A"| . 

v. u y — > tana as y — > oo and hence \u y \ is bounded. 

Since u yy (0) > 0,[ii]shows that u y is either strictly monotone increasing or has 
one strict local maximum. Therefore Irnl and [rvl yield the claim. 

vi. u yy has at most one zero. Therefore u yy — > as y — > oo. 

This follows from [n] and [v] 

To prove theorem [I] we show first that the asymptotic angle lim y _ s . 00 a (y) of 
u goes to when C goes to or oo in ©. Together with continuous dependence 
of the solutions u on C > this is enough to prove theorem [TJ 

As a first step we show that the first derivative of the solutions u blows up, 
when C goes to or oo in ©. In view of |v] this is useful for investigating the 
behaviour of the corresponding asymptotic angles. 

Lemma 2.2. The solutions u of (0|) satisfy supj — > oo if C —¥ or 

C -> oo. 

Proof. First we treat the case C — > oo. We can write © as 

n — 2 u yu y 

+ 2 " T" + l + (Uy) 2 ' 

Since u(0) = C, [n] implies that u (j/q) — > oo as C — > oo for any fixed j/o > 0. 
Therefore the last equation yields either u y (j/o) — > oo as C — > oo or u yy (j/o) — > 
oo as C — > oo. In both cases we can deduce the claim. 

For the case C — > suppose not, so we may assume without loss of generality 
that sup[ 0iOO ) u y < D, D > for all C > small enough. First we use the 
clearing out lemma (see appendix to show that u must be large compared 
to such C > away from 0. Therefore assume that for C > small enough 

u (LC) < KC, K > 0, L > 

for a constant K, where L will be determined later. 

Now we want to apply the clearing out lemma (see appendix [2) to the self- 
expander M (u) (considering u as a function on M here) with respect to the ball 
Blc (0, C, 0, . . . , 0) C R n . Estimating the Hausdorff measure yields 

•H™- 1 (M (u) n B LC (0, C, 0, . . . , 0)) < DK n - 2 LC n - 1 

1 

where now L> = D{n). Now choose L (D, K, n) > 0, L = K 2 , to make 

DK n ^ 2 LC n ^ 1 < eoi™ _1 C" _1 , where £o is the constant from the clearing out 
lemma. 

Therefore the clearing out lemma says that the clearing out time tc of M(u) 
with respect to B L q/4 (0, C, 0, . . . , 0) can be estimated as tc < cL 2 C 2 for some 
constant c > 0. But since M[u) moves self-expanding we also have ^/T+TcC > 
C + I£ and hence tc > 4- - a contradiction. 
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This means u (LC) must be large for any C > small enough. More precisely 
there must be a sequence (Cfe) feeN which satisfies 



lim LC k =0, u k (LC k ) > kC k . 

k— too 



where Uk are the associated solutions of ([5]) with Ufe(0) = Ck and L(D,k,n) as 
above. By [n] we have 

rLCk 

u k (LC k ) =C k + / {u k )y dy <C k + LC k max (u k ) 

JO [0,LC k ] 

and therefore max [0jic j (u k ) y > jj, D = D (n) > 0. 

Now we use continuous dependence to get contradiction to the last estimate. 
For any fco € N we can choose yk > sufficiently small to make (uk ) y < ^ 
on [0, j/fc ]. Since u y depends continuously on the initial values of © we can 

choose fco e N such that maxr ^ i (uk ) y - {uk) v < 



we get 



(Uk)y < {Uk a )y-(Uk) y + (U k0 )y 



for k > fco- Therefore 
1 



< 



D 



on [0,yfc o ] for any k > fco - a contradiction. 



□ 



From the last lemma we know that u y blows up somewhere on (0, oo) as 
C — > or C — > oo. To get the desired behaviour of the asymptotic angles, it is 
therefore enough to show that u y does not decrease too much after a possible 
maximum. 

Lemma 2.3. For solutions u of 0) we have a — > ^ as C — > or C — > oo. 

Proof. Because of[vi]and lemma [2?2l we may assume that for any C > there is 
precisely one zero of u yy , yc > 0. We have from © for y > 



u 

y 



2(n-2) 



2it,, 



uy 



y (l + K) 2 ) 



y 2/ 



i + Kr 



This means (-) > on [yc, oo) for any C > 0. By[m]we have - — »• 



as y — > oo, so it is enough to show ^ — >■ oo at some point in [yc, oo) if either 
C -> oo or C -> 0. 

We treat the case C — > oo first. Suppose not, so without loss of generality 
there exists D > such that "^f' 1 < 13 for any C > large enough. Hence by 

the above equation — u(^ c ^c Uy — ^ an< ^ lemma 12721 (yc) — > oo as 
C — > oo. Therefore we have u (yc) yc — ► as C — ► oo. But by[n]we have u> C 



on [0,oo), which then yields 



"(j/c) 
yc 



> C 2 - a contradiction. 
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For the case C — > again suppose not. Without loss of generality we have 
then | < /) on [yc, oo) for some D > and for any C > small enough. 

For such C > we get < £> an d 

u (2yc) = u (yc) + / Uydy > y c u y (2y c ) ■ 
Jyc 

Hence u y (2yc) < 21?. On the other hand © written as above in the proof 
yields 

u (ggg) < 2(n - 2) u, (2gg) 
2y c " 2y c u(2y c ) 2 

Therefore we get u (2yc) yc > j) for C > small enough, after possibly ad- 
justing D > 0. So it remains to show yc — > as C — > in order to get a 
contradiction. 

Assume that yc as C — > 0, so without loss of generality there exists 
e > such that yc > £ for any C > small enough. As in the case C — > oo we 
get from u v (yc) — > oo that ycu (yc) — > as C — > 0. 

In view of the bound on "^f ^ and[ii]we see that w — ► uniformly on [0,e] 
as C — > and 

vc\ , /"^ j~ ^.vc ( yc 



u (yc) = u l^— j + u y dy > —u y ^ o 

Therefore u a f 1 < ^ an d especially u y < D on [0, |] , after possibly adjusting 
D > 0. 

We can now apply Brakke's clearing out lemma (see appendix) for M(u) 
(again with u : K — » M) with respect to the ball Be (0, C, 0, . . . , 0) to get a 
contradiction. Estimating the Hausdorff measure yields, using our bound on u y 



1 (M(u)nBe (0, C, 0, . . . , 0)) <D(n)ema,xu 



n-2 



Since max|- £i u n 2 — >• as C — > the clearing out time ic of M(u) with respect 

to Be (0, C, 0, . . . , 0) therefore satisfies < c x for C > small enough, where 
c > is a constant. But M(u) moves self-expanding, so y/tc~TTC > C + § 

~2 

which means tc > g | g ^ - a contradiction. 

□ 

Finally the next lemma is the desired stability result for ([6]). 

Lemma 2.4. The asymptotic angle of solutions of {6p, a — a(C) : (0, oo) — > 
(0, §), is a continuous function of the initial condition C > 0. 

Proof. Assume not. Then there exists e > 0, Co > and a bounded sequence 
(Cfc) fcgN in (0, oo) with Ck — > Cq as k — > oo and jtanac 1 ,. — tanacj > £ for all 
k e N. The associated solutions of ^ are denoted by (wfc) fe6N and uq. 
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From Irvl we see that < sup y>0 (uk) y < D for all k e N, D > a constant. 
In the following we will adjust the constant D > implicitely as necessary. In 
view of[vI]we assume that (uo) yy has precisely one zero yg. The case (uq) > 
on [0, oo) can be handled in the same way. 

Since we have continuous dependence on the initial conditions on compact 
intervals we may assume that all Ck are close enough to Co, so that each (uk) yy 
has precisely one zero close to yo- Therefore we get sup y>0 {uk) yy < D for all 
k £ N. Differentiating ^ with respect to y > yields 



— (n — 2)u y 



yUyy 



1 



i + (u y y 



(3) _ (Uyy) 



i + (u y y 



In view of this means (uk) vv can not have arbitrarily small local extrema. 



Therefore 



< D for all k £ N. 



sup y > (u k ) yy 
Hence we can estimate for y > 0, using HI and IvTl 

'uk\ 1 

y 



y 



V 

' n- 2 



y 



(«fc). 



y 



(Uk) 



mi 



i + (u y y 



< 



D 



y 



2 ' 



So by |m] for any S > there exists ys > with 



y 



tan q;^ 



< <5, 2/ > ys, k e N. 
Using our assumption on aCo, ac t we can therefore choose 8 > such that 



> 



We can write the difference of the ODEs for uo and Uk at y > as £i 
6+6+6, where 



6 = - («o - Wfc) 

y 

6 = - (U k )y 

Myy ( U k) yy \ 



y\i + (u )l i + K) 2 



y , 

2 /n-2 n-2^ 



y \ u k uo J 

By the previous considerations we can find y > y$ such that £3 < | and 
6 < § for any k £ N. By continuous dependence we can find fc G N such that 
6 < §• Therefore we get |6 + 6 + 61 < § - a contradiction. 

□ 
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Putting the last three lemmas together we arrive at theoremQ] a* rit £ (0, §) 
follows from [v] and from the clearing out lemma (see appendix [S] and |AIj ) . As 
stated in lemma 12.11 the existence of a two-sheeted solution (for any cone angle 
a e (0, ■§)) is asserted using results of Ecker and Huisken (see [EH]) which 
ensure existence of self-expanding evolutions of the two Lipschitz graphs C a 
and C-ji — q. 

Additionally standard theory of differential inequalities shows that for C > 
y/2 the asymptotic angle a (C) is strictly monotone increasing for solutions of 
©. We believe it is everywhere strictly monotone, apart from C cr u, a (C cr i t ) — 
a* crit . This would imply that for each a* rit < a < | there are precisely two 
self-expanding, smooth evolutions of D a of the form M(u). In general there 
might be more self-expanding evolutions of D a , possibly also non-rotationally 
symmetric (see |AIj ) . 

Remark 2.5. Solutions of ^ are stationary for the functional 

K[M} = J^p^pjdH^ix). 

The authors of [All sketch a proof for the existence of one- sheeted self-expanders, 
asymptotic to D a which are minimizers o/K. Furthermore they indicate how 
one might prove a version of theorem Q] using this approach. 



3 Touching fluid droplets 

As mentioned before we want to apply the previous results to study the be- 
haviour of touching fluid droplets. These are assumed to have locally conical, 
rotationally symmetric shape, i.e. the shape of a smoothing of D a . The follow- 
ing definition makes this formulation precise. From now on we set the dimension 
to n = 3. 

Definition 3.1. We call M a — M (u a ) a smoothing of the double cone D a with 
angle < a < ^ , 7 = tana, if for a < < b 



u a (y) 



l\y\ ify<aorb<y 
s{y) if a<y <b 



such that s > and u a e C 2 '^ (K) , [3 > 0. M a is said to lie outside of the 
double cone D a if s (y) > j\y\ for a < y < b. u is the generating function of 

Ma- 

Clearly any smoothing of a double cone stays rotationally symmetric under 
mean curvature flow. One can compute (see [SI]) that ([1} for a mean curvature 
flow evolution M t of M a is equivalent to 



9 

at , , / » \" u 



8 ih~ u 1 

° u = ^ s - - (7) 



(*«)■ 
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where u — u(-, t), u(-, 0) = u a generates M t . 

For any smoothing M a we have short-time existence of a solution M t of |T]) 
on a maximal time interval [0,T), T > 0. Furthermore the solution must be 
smooth for t > and every finite time singularity must be due to pinching, i.e. 
infgu(-,i) — > as t — > T. This holds even without any growth assumption on 
the initial generating function (see [Sj and |LSj). 

In fact a sphere comparison argument (see [EC] ) shows that 

lim mxf\u(-, t) — > as t — > T 

t^T R 

must hold for finite time singularities. 

This agrees with intuition about repulsion (pinching in finite time) and co- 
alescence (long-time existence) of fluid droplets. 

Using techniques by Eckcr and Huisken (see |EH| ) one can derive global 
height estimates for rotationally symmetric solutions of the mean curvature flow 
which yield, using the results from [BBj . the following comparison principle (see 
[BO]). 

Lemma 3.2. Let M ai , M a2 be two smoothings of the double cone with u ai < 
u a2 for the associated generating functions. Denoting the generating functions 
of the two evolutions with u 1 and u 2 we have then u 1 < u 2 as long as the 
solutions exist. 

Corollary 3.3. Any smoothing of the double cone M Q , < a < \ has a unique 
evolution by mean curvature. 

Using these results we can now define what coalescence and repulsion mean 
within our model. 

Definition 3.4. An angle < a < ^ is called a repulsion angle if there exists 
a smoothing of the double cone M a which is outside of the double cone and such 
that the mean curvature flow evolution of M a pinches in finite time. 

Lemma 3.5. An angle < a < ^ is not a repulsion angle if and only if there 
is a smoothing of the cone M a , for which the evolution under mean curvature 
flow exists for all t > 0. 

Proof. By definition any angle < a < ^, that is not a repulsion angle, must 
have a smoothing M a for which the evolution under mean curvature flow M t 
exists for all t > 0. 

So suppose for < a < 5 there exists a smoothing M Q , such that the 
evolution M t exists for all t > 0. Let M a be another smoothing, that is outside 
of the double cone. We denote its evolution by M t , t € [0,T). 

We know that the mean curvature flow is invariant under parabolic rescaling 

x i y Xx, 1 1 — ^ X 2 t 

for any scaling parameter A > 0, x € M t and t 6 /. Let M t A be the rescaling of 
Mt- Note here that any double cone D a is invariant under the scaling x i— > Xx. 
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Since M a is outside of the double cone we can therefore choose A > sufficiently 
small in order to get initially 

< u a 

for the corresponding generating functions. By lemma 15721 Mt must exist for all 
t > 0, therefore a is not a repulsion angle. □ 

In view of the last lemma we make the following definition. 

Definition 3.6. An angle < a < \ is called a coalescence angle if it is not a 
repulsion angle in the sense of definition \3.4\ 

Given a field strength (respectively a cone angle) at which coalescence oc- 
curs any lower field strength (respectively greater cone angle) should lead to 
coalescence as well. The same should hold for repulsion angles with higher field 
strength. The next two lemmas shows that this is true within our model. 

Lemma 3.7. Let ao be a coalescence angle. Then for any angle a > ao and 
for any smoothing M a that is outside of the double cone the mean curvature 
flow with initial data M a exists for all t > 0. Therefore any angle a > ao is a 
coalescence angle. 

Proof. As in the proof of lemma 13.51 we can scale appropriately and then use 
the comparison principle from lemma 13.21 □ 

Lemma 3.8. Let ao be a repulsion angle. Then for a < ao any smoothing of 
the double cone M a with angle a must pinch in finite time. This means any 
angle a < ao must be a repulsion angle. 

Proof. Again, as for lemma 13.51 we can scale appropriately and then use lemma 
I3~21 □ 

The observations from the experiments in RB suggest the existence of a 
critical angle for the behaviour of the system. As the next lemma shows, this is 
also true for our model. 

Lemma 3.9. There is a critical angle < a cr u < § such that any angle 
a < cicrit is a repulsion angle and any angle a > a cr it is a coalescence angle. 

Proof. Let 



a = sup |0 < a < — : a is a repulsion angle | 

a = inf|o<a;< — : a is a coalescence angle | . 

According to lemma [3751 any cone angle a < a is not a coalescence angle. This 
shows that a < a. Clearly any angle is either a repulsion or a coalescence angle, 
hence we have a = a = a cr u- D 

Using the constructed one-sheeted self-expanders we can now determine 
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3.1 Determining a cr u 



Lemma 3.10. Any double cone smoothing with angle a > a* rit (3) that lies 
outside of the double cone has an evolution which exists for all t > 0. 

Proof. Let M a be a smoothing of D a with a > a* crit {3) and generating function 
u a . Let s be the generating function of the self-expander asymptotic to D a * ( 3 ) 
which exists by theorem [TJ As in the proof of lemma [3~31 we can do a parabolic 
rescaling and therefore assume s < u a . Then the evolution of M a must exist 
for all t > by the comparison principle, lemma 13.21 since every finite time 
singularity must be due to pinching. □ 

Lemma 3.11. Any double cone smoothing with angle a < a* rit (3) must pinch 
in finite time. 

Proof. We follow a level-set flow argument from |AIj . 

Assume that there is a double cone smoothing M a , a < a* rit (3), such that 
its evolution by mean curvature M t exists for all t > 0. 

Let T t , t > 0, be the level-set flow of D a . T t can be characterized as follows. 
M 3 \rt is the union of all level-set flows At such that Ao is compact and lies in 
R 3 \D a . 

First wc show that € L. So let At be a level-set flow with Ao C M. 3 \D a 
compact. Ao must either lie in the convex hull of one of C a or C v - a or outside 
of D a . Using the maximum principle for one level-set flow and one smooth flow, 
we see that in the first case Ao is pushed away from by the graphical self- 
expanders of Ecker and Huisken (part (i) of lemma l2~Tj) and therefore A t . 
In the second case we can parabolically rescale M a as in the proof of lemma l3~5l 
and therefore assume Ao f~l M a = 0. Then we can apply the maximum principle 
to see that At. Hence by the above characterization we must have G r t . 

Ti is rotationally symmetric, so Ti = M(X) for some closed set X C K 2 . 
In fact the boundary dTi is smooth and d X consists precisely of curves of the 
type in lemma |2~T1 (see jAIj). 

Since € T\ there must be a curve of type (ii) in lemma I2TT1 in dX, which 
is asymptotic to a a U a n - a . This yields a smooth, rotationally symmetric, 
one-sheeted self-expanding evolution of D a - a contradiction. □ 

In view of the definition of a cr u the last two lemmas yield the following 

Corollary 3.12. a cnt = a*„ t (3). 

3.2 Conclusions 

Using a model based on mean curvature flow, we obtain a critical cone angle 
a crit for the behaviour of oppositely charged droplets of fluid. More precisely 
this means any smoothing of the double cone (see definition 13. II) with angle less 
than a cr u must pinch in finite time (repulsion). Assuming a smoothing has an 
angle greater than a cr n and lies outside of the double cone, its evolution must 
exist for alH > (coalescence). 
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We can show that ct cr n — ct* rit (3), where ct* rit (3) is the critical angle for 
the existence of smooth, rotationally symmetric, self-expanding, one-sheeted 
evolutions of double cones and a* crit (2>) ss 66°. This coincides with observations 
from experiments (see |RB] ) which predict a critical angle of 60° — 70°. 

Finally we want to compare our model with the one in |RBlj . For that 
approach it is assumed that the bridge between the touching droplets minimizes 
area under a volume constraint. This corresponds to constant mean curvature 
surfaces of revolution (Delaunay surfaces) which are fitted to linear double cones, 
similarly as in definition 13.11 (unlike in definition 13.11 the associated generating 
function is only continuous). The associated capillary pressure p is assumed to 
determine the behaviour. 

The critical shape has p — which yields a rescaled catenoid. Suitable 
smoothings of the associated surfaces with p > provide double cone smoothings 
in the sense of definition 13.11 which pinch in finite time. This is in agreement 
with RBI . However lemma 13.111 shows that down to a certain po < the 
(smoothings of the) associated surfaces with p < still pinch in finite time 
under mean curvature flow. Therefore our predicted critical angle is greater 
than the predicted critical angle from [RBI] , which is approximately 59°. 

A The clearing out lemma 

This version of the clearing out lemma is due to Brakke. It says that if the 
area ratio of a piece of a solution of mean curvature flow in a ball is initially 
small, then it must leave the ball with one quarter of the radius after a time 
proportional to the initial area ratio times the radius squared. Proofs can be 
found in [BR] or [EC] . 

Theorem. Let (M t ) t>0 be complete, properly embedded hyper surfaces in 1", 
n > 2, evolving by mean curvature flow. There exist constants c(n) > and 
£o(n) > such that if Mq satisfies 



K 1 - 1 (M n B p (so)) < e p 



,n-l 



for some xq £ W l and p > 0, then 




for t = ce p 2 . 
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